Thermal fluctuations in macroscopic quantum memory 
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We describe macroscopic quantum memory devices based on type-II toroidal superconductors and 
estimate in one case and compute in another the rates at which quantum information stored in 
these devices "degrades" because of thermal fluctuations. In the case when the entire solid torus is 
superconducting, the Boltzmann factor in the rate corresponds to a well-defined critical fluctuation, 
and the rate is suppressed exponentially with the linear size of the system. In the case when 
superconductivity is confined to the surface of the torus, the rate is determined by diffusive motion 
of vortices around the torus and does not depend exponentially on the linear size; we find, however, 
that when the two dimensions of the torus are comparable the rate does not contain the usual 
volume enhancement factor, i.e. it does not grow with the total surface area of the sample. We 
describe a possible way to write to and read from this quantum memory. 
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Quantum memory is a device capable of reliably stor- 
ing linear superpositions of quantum states. It will be 
a part of quantum computer when (if) that latter is fi- 
naUy buih and may be useful for other applications as 
well. (For a recent review of quantum computing with 
an emphasis on fault tolerance see ref. ||l|.) 

To work as quantum memory, a physical system must 
satisfy a number of requirements. First, it must have 
at least two fairly stable quantum states. These states 
form a basis for linear combinations that can be stored 
in the device. For example, the basis may be formed 
by perturbative quantum states built near local energy 
minima, and stability of the basis states may be ensured 
by a large potential barrier separating them. In such 
cases we will loosely refer to the basis states as the ground 
states, or vacua, even though these ground states may not 
be degenerate in energy and in some cases may contain 
localized excitations. We note, though, that for some 
purposes it may be desirable to have ground states that 
actually are (nearly) degenerate in energy. If two basis 
states \ipi) and IV'2} forming a linear combination 
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are degenerate, the ratio C2/C1 will be preserved by the 
evolution. When the basis states are not degenerate, the 
relative magnitude of ci and C2 will be preserved, but not 
the relative phase. The relative magnitude, however, can 
be arbitrary. In comparison, a classical two-state system 
will only store two values, referred to as and 1. 

For long-time quantum storage, one will probably need 
to build in some redundancy, so that the basis states refer 
to many microscopic (local) degrees of freedom. However, 
redundancy is helpful in protecting quantum information 
only when the local degrees of freedom in the basis states 
are sufficiently entangled, i.e. the basis states cannot be 



identified by local measurements. This condition rules 
out, in particular, any system in which a ground state de- 
generacy is due merely to spontaneous symmetry break- 
ing by a local order parameter. (To see why, consider an 
easy- axis magnet, in which magnetization can be in one 
of two directions. The direction of magnetization can be 
found by measuring local magnetization in a relatively 
small region.) The reason why entanglement of local de- 
grees of freedom is necessary for long-time quantum stor- 
age is that local measurements will in effect be performed 
by external noise, and if they can indeed distinguish be- 
tween the basis states they will destroy the stored quan- 
tum information (the Schrodinger-cat scenario). 

In a real sample, tunneling transitions between the 
basis states will cause quantum memory to deteriorate. 
Nevertheless, if the basis states are sufficiently entan- 
gled, tunneling between them will have to involve many 
local degrees of freedom, and the tunneling probability 
will be strongly suppressed. Known examples in- 
clude fractional quantum Hall and similar types of rigid 
ground states on tori. In these typical tunneling 

fluctuation consists of creating a vortex-antivortex pair, 
transporting the vortex and the antivortex around the 
torus, along topologically distinct paths, and then anni- 
hilating the pair. It has been argued that, at zero temper- 
ature, the tunneling probability (and the associated en- 
ergy splitting between the ground states) is generically of 
order exp{—L/l) where L is the size of the system, and 
I is some correlation length. So, the zero-temperature 
tunneling should not be a problem in practice, as long as 
one can keep the size of the system sufficiently large. 

Of more concern are thermal fluctuations. At finite 
temperature, there will be a sea of vortex-antivortex 
pairs, with density proportional to exp(—Fo/T), where 
Fq is the free energy of a single vortex. (This assumes 
that the temperature is still low enough, so that no phase 
transition occurs.) One expects that motion of these 
"preexisting" vortices can effect transitions between en- 
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sembles built near different ground states. The question 
is, then, what is the rate of such transitions, at a given 
temperature T. A single transition is sufficient to de- 
stroy the stored quantum information. So, the rate of 
the transitions will also be the rate at which the quan- 
tum information "degrades" , and the corresponding time 
will be an estimate for the maximal duration of reliable 
storage. 

In this paper, we will compute in one case and estimate 
in another the rates of finite-temperature transitions be- 
tween ensembles corresponding to different ground states 
for some of the simplest systems exhibiting multiple 
ground states and macroscopic entanglement. One sys- 
tem is a type-II superconducting film grown on the sur- 
face of a torus. In Sect. 2 we review the origin of multi- 
ple ground states in a type-II superconductor on a torus. 
The presence of multiple ground states in this case can 
be seen either via manipulations with vortices and single 
electrons, which produce a nontrivial phase when trans- 
ported around each other, or via a semiclassical argu- 
ment. Transitions between different classical vacua are 
topological transitions, which change a winding number 
of the gauge and Higgs fields. In Sect. 3 we construct 
a correlator that measures the rate of topological transi- 
tions at finite temperature. This correlator is analogous 
to the one proposed in ref . jj] to measure the rate of topo- 
logical transitions in the electroweak theory. In Sect. 4 
we compute the rate. The main ingredient of the com- 
putation is that vortices are well separated, and their 
motion is diffusive, i.e. associated with a large viscosity. 

Our main result, for a film of fixed thickness, is that 
although the rate of topological transitions is indeed pro- 
portional to the vortex density, and so is not suppressed 
exponentially with L (the size of the system), there is 
a power-law suppression. This suppression can be de- 
scribed by saying that there is no volume enhancement 
of the rate, i.e. as long as the two dimensions of the 
torus stay comparable, the rate will not grow with the 
total volume (while the total number of vortices of course 
will). Equivalently, the rate per unit volume will decrease 
with the total volume. This absence of macroscopic en- 
hancement is directly related to the diffusive nature of 
the vortex motion. 

It is easy to redesign the device so that the sup- 
pression of the finite-temperature rate becomes expo- 
nential with L. Imagine making the superconducting 
film thicker, so that vortices resolve into Abrikosov flux 
lines; the free energy of those grows linearly with their 
length. In the limiting case, which is the second sys- 
tem we consider, the entire solid torus is superconduct- 
ing, and a topological transition is mediated by a well 
defined critical fluctuation — a critical flux line, whose 
energy is proportional to L. The Boltzmann factor in 
the rate is exp(— i?o/T), where now Eq cx L, so the 
finite-temperature rate is suppressed exponentially with 
L. The zero-temperature tunneling rate is suppressed 
even stronger, as an exponential of L^. So, a solid su- 
perconducting torus (or a wire, or a ring, or a hollow 



cylinder) is a good candidate for stable quantum mem- 
ory. In the concluding section we discuss a possible way 
of writing quantum information to and reading it from 
this device. 

We nevertheless retain interest in the two-dimensional 
case (the film), because a universal quantum computa- 
tion is theoretically possible with non-Abelian anyons [pi , 
and systems in which those have been argued to occur ||5[| 
are two-dimensional. In the concluding section we also 
discuss whether our results teach us anything about these 
more complex cases. 

II. GROUND STATES OF A TOROIDAL 
SUPERCONDUCTOR 

Existence of multiple ground states in a type-II super- 
conductor on a torus can be deduced from the presence 
of two types of local excitations, vortices and single elec- 
trons, with their corresponding values of flux and charge. 
It can also be obtained from an explicit semiclassical con- 
struction of the ground states. In this section, we use 
the Ginzburg-Landau (GL) theory for description of the 
ground states. We interpret the GL expression for energy 
as an effective Hamiltonian for slow degrees of freedom 
(rather than as a thermodynamic potential, like free en- 
ergy). So, we treat the GL fields as quantum fields. 

The GL Hamiltonian of a superconductor is 

H = J d^x (C|(V + igA)i;\^ - a|V'P + b\ip\^) + Hem , 

(2) 

where C, a, and b are positive coefficients, 

9 = 2e/c , (3) 

2e is minus the electric charge of a Cooper pair (e > 0), 
and c is the speed of light; h = 1 everywhere. We con- 
centrate on the extreme type-II case; the corresponding 
condition on the parameters is 

« b . (4) 

The Hamiltonian of electromagnetic field is taken, for 
simplicity, in the relativistic form: 

i^EM = ^ Jd^xiE^ + B^) . (5) 

In (^, (H) "0 is the complex "order parameter" field (it 
is not really an order parameter because it is not gauge- 
invariant [|6|), A is the electromagnetic vector potential, 
E and B are the electric and magnetic fields. 

We first consider a superconducting film that extends 
from z = to z = d in the z direction and is peri- 
odic (toric) in the x and y directions. These periodic 
boundary conditions define what may be called a "math- 
ematical" torus, as distinct from the surface of a physical 
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torus, a "doughnut" , that one may produce in a labora- 
tory. Later, we will discuss the distinction in more detail 
and will also consider the case when the entire solid torus 
is superconducting. 

The vortex of the theory (^) is a short (of length d) 
Abrikosov flux line whose axis is parallel to the z axis. 
A vortex carries magnetic flux of 27r/g — irc/e. So, if we 
break a Cooper pair and transport one of the electrons 
around the vortex, the wave function of the system will 
acquire a nontrivial Aharonov-Bohm factor of —1. As 
shown in refs. 0, whenever transport of local excita- 
tions around each other produces such a nontrivial factor, 
the ground state of the system on a torus is degenerate, 
up to an energy splitting decreasing exponentially with 
the system's linear size. For the present case, it comes 
out that the ground state degeneracy on a torus is at least 
four-fold. We do not reproduce the argument here, as it 
can be found in the above papers. Besides, in our case 
the vacuum structure admits a semiclassical interpreta- 
tion, which allows us to obtain all the requisite results in 
a different way. 

Consider classical vacua of (||), i.e. configurations of 
the lowest energy. On a torus, these are: 

2tt ( n^e^ n,,e,.\ 
A - — ^ + , (6) 

t/j ^ t/joexp {-2TTinxx/Lx - 2TTmyy/Ly) , (7) 

where and Uy are arbitrary integers, and Sy are 
unit vectors in the two directions, and L^, Ly are the 
corresponding dimensions of the torus; 

% = (a/26)i/2 . (8) 

The integers and ny are the winding numbers of the 
configuration: they count how many times the phase of 
-0 winds as one travels along the torus's noncontractible 
loops. We consider the case when and Ly are compa- 
rable and assume, for definiteness, that 

Lx> Ly , (9) 

i.e. that the larger loop of the torus is in the x direction. 

Tunneling processes mix the perturbative vacua built 
near the configurations (|^)~(0) into linear combinations, 
6*- vacua, analogous to those of the four-dimensional QCD 
1^. If we denote the perturbative vacua as \nx,ny), the 
0- vacua are 

\^x,Oy) = ^ exp(i6' ny) , (10) 

where 9x and 9y run from to 2tt. In this case we need 
two 9 angles because there are two winding numbers, Ux 
and Uy. A more important difference from QCD, though, 
is that in the present case the tunneling amplitudes, and 
hence the energy splittings among the 9 vacua, are ex- 
ponentially suppressed with Lx or Ly. This exponential 
suppression was found in ref. in a slightly different 



context, see also ref. [^. It can be explained as follows. 
A typical tunneling fluctuation consists of a vortex and 
an antivortex, which travel along topologically distinct 
routes: the vortex travels distance AL, and the antivor- 
tex distance Ly — AL (if wc consider transitions that 
change Hx). At least one of these distances is macro- 
scopically large, and to travel that far the object has to 
move very fast, or to stay in existence for very long, or to 
achieve a good balance between these two extremes. One 
finds |] that even the fluctuation that achieves the opti- 
mal balance still has a Euclidean action proportional to 
Ly, resulting in an exponentially suppressed amplitude. 

In what follows we will assume that system is suffi- 
ciently large, so that the tunneling processes that change 
Ux and Uy are practically nonexistent. In this case, the 
linear combinations ( p^ are no longer special, and an 
equally good basis in the ground state subspace is pro- 
vided by the perturbative vacua \nx,ny) built near the 
classical solutions (|6|)-(0). From the nontrivial proper- 
ties of excitations, with respect to transport around each 
other, we have learned that, when tunneling is neglected, 
there are at least four degenerate ground states. Now we 
find infinitely many degenerate vacua |na;,ny). It is easy 
to make four from infinitely many. Note that the ground 
states \nx,ny) and \nx + i,ny) can be distinguished by 
breaking a Cooper pair and transporting one of the elec- 
trons around the torus in the x direction. Say, for Ux = 
the electron will pick no phase factor, while for rt^; = 1 it 
will pick a factor of —1. On the other hand, given that 
the charge of electron is the minimal charge in the sys- 
tem, there is no way to distinguish between \nx,ny) and 
\nx + 2,ny). Similarly, one cannot distinguish between 
\nx,ny) and \nx,ny + 2). So, in the absence of tunnel- 
ing, instead of the infinitely many vacua \nx^ riy) we may 
as well consider only four "equivalence classes" , corre- 
sponding to Ux and Uy both being even, one being even, 
the other odd, and both being odd, respectively. The 
four vacua deduced from the quantum numbers of the 
excitations are representatives of these four equivalence 
classes. 

Although, as we have seen, in the absence of tunneling 
we do not have to consider the entire infinite "lattice" of 
the vacua \nx,ny), sometimes it is convenient to do so. 
In particular, in the next section we will see that ther- 
mal fluctuations in the winding numbers are conveniently 
viewed as diffusion of Ux and Uy over an infinite lattice 
made by pairs of integers. 

Now consider a type-II film that sits on the surface 
of a solid torus, a "doughnut", whose bulk is not su- 
perconducting. There are still two winding numbers, nx 
and ny. For example, Hx in this case is simply the to- 
tal magnetic flux through the doughnut's hole, in units 
of the flux quantum $o = 27r/(7. One can change Ux 
to [ux -l- 1) by dragging an extra flux quantum from the 
outside, through the bulk of the doughnut. This is equiv- 
alent to creating a vortex and an antivortex on the outer 
side of the doughnut, transporting them along topolog- 
ically distinct paths to the inner side, and annihilating 
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them there, cf. ref. ||^. In quantum theory, this pro- 
cess occurs spontaneously, as a quantum fluctuation. It 
is a tunneling process between two distinct ground states 
that differ by one unit of rix- The conclusion that the 
tunneling rate is suppressed exponentially with Ly (for 
transitions that change rix) still applies. 

One can switch between the ground states "by hand" , 
i.e. by dragging appropriate fluxes with the help of exter- 
nal solenoids. Switching from jV^i) to \tp2), for a system 
that was initially in the linear superposition is equiv- 
alent to interchanging ci and C2. It is hard to say, though, 
if this "quantum switch" can serve any useful practical 
purpose. 

On a "doughnut" , the superconducting current in a 
ground state with n^^ ~ 1 and Uy = is of order 
c^q/C, where C is the self-inductance of the device: 
C ^ Lxln{Lx/ Ly); in what follows we assume the log- 
arithm here to be of order one. The current density is 
then inversely proportional to the surface area of the sam- 
ple. Because the current density is so low, ground states 
differing by a few units of Ux or Uy are practically indis- 
tinguishable with regard to energies of local excitations. 
This circumstance plays an important role in preserva- 
tion of quantum coherence at finite temperature. 

On the doughnut, as opposed to the "mathematical" 
torus — a rectangle with periodic boundary conditions, 
the ground states corresponding to different values of 
and Hy are not exactly degenerate even in classical the- 
ory: for different values of Ux there are different amounts 
of energy associated with the magnetic field trapped in 
the doughnut's hole. This energy has very little influence 
on the rate of topological transitions, so calculation of the 
rate can be carried out on the "mathematical" torus. On 
the other hand, a real device will be a "doughnut" , and 
in that case the magnetic energy will lead to discrete (la- 
belled by Ux) energy levels. A resonator tuned to the 
energy difference between two such levels may then be 
able to write linear superpositions of quantum states to 
this device, or to a solid superconducting torus, which we 
discuss later. Estimates related to this writing technique 
are given in the concluding section. 

III. TOPOLOGICAL TRANSITIONS AT FINITE 
TEMPERATURE 

Consider, for definiteness, two ground states that differ 
by one unit of Ux and have the same Uy. These states 
may play the role of the basis states {ipi) and \ip2) forming 
the linear combination (|l|). When the system is at a 
finite temperature, excited states are also occupied. For 
low-lying, perturbative excited states, we can distinguish 
between states built near lipi), let us call them IV"!"'): 

and those built near |'02), call them \ip^2^''); we neglect 
tunneling between these two sectors. It is essential that, 
for superconducting tori that we consider here, excitation 
energies are nearly the same for the two sectors, even 



when (as on a "doughnut") there is a sizable difference 

(n) 

between the corresponding ground state energies. If -E^ 

denotes the energy of 1^1"^}, E^°^ the energy of [i — 
1,2), and the excitation energies are, to a good accuracy, 
equal: 

there is a quasiequilibrium state described by the follow- 
ing density matrix: 

p-5]«;„|V'("))(V'(")l, (12) 

n 

where 

M")) = ci|^("))+C2|4"^) , (13) 

and w„ = 7Vexp(-A£;(")/r); M is a normalization fac- 
tor, and T is temperature. This state is no less quantum- 
coherent than the ground state (p. 

The density matrix (|l^) is not a fully equilibrium one, 
because it does not include occupation of higher excited 
states, those that cannot be regarded as being "near" 
either \'tpi) or |V'2)- Thermal occupation of these states 
will determine the rate of transitions between the two 
sectors. Our goal in this section is to set up formalism 
for calculating these transition rates. We begin with the 
case of a type-II film on the surface of a "doughnut" . 

Any two configurations that differ by one unit of Ux 
or one unit of Uy are separated by a potential barrier 
whose height is, to a good accuracy, twice the energy of 
a static vortex. In general, a system at finite temperature 
does not need to tunnel under a barrier; it can go over 
it as a result of a thermal fluctuation. In many cases, 
the rate of these thermal transitions can be computed by 
considering vicinity of the fluctuation corresponding to 
the top of the barrier This fluctuation is called the 
critical fluctuation. For a toroidal type-II film, however, 
calculational schemes based on expanding near a critical 
fluctuation are completely useless, for the following rea- 
son. The top of the barrier in this case corresponds to 
a vortex and an antivortex separated by distance Ly/2 
(for transitions that change Ux), see Fig. 1. But at a 
finite temperature there is a finite density of vortices and 
antivortices, with a typical distance between them that 
is much smaller than Ly/2. In this situation, a pair of 
widely separated vortex and antivortex cannot have any 
special significance. Accordingly, we expect that the rate 
of topological transitions will be determined by motion 
of vortices already populating the medium. Nucleation 
and annihilation of vortex-antivortex pairs will merely 
maintain the equilibrium concentrations of vortices and 
antivortices. 
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FIG. 1. Critical fluctuation, representing the top of the po- 
tential barrier between neighboring classical vacua, in the case 
when only the surface of torus is superconducting (shaded 
area). The critical fluctuation consists of a vortex and an an- 
tivortex lying on a cross-sectional diameter. It "melts" at a 
finite temperature, because many vortices intervene between 
these two, and so plays no role in thermal transitions between 
the vacua. 

In contrast to the two-dimensional case (film), a criti- 
cal fluctuation can be readily identified at finite temper- 
ature in a solid superconducting torus, or a loop of thick 
superconducting wire. A solid type-II torus has multiple 
ground states, although not as many of them as a torus in 
which superconductivity is confined to the surface. Loops 
in the y direction are contractible through superconduc- 
tor, so there is no winding number Uy that would corre- 
spond to those. But Ux still exists and still counts the 
number of flux quanta trapped inside the loop. Chang- 
ing Ux by dragging a flux through the loop is still oper- 
ational, but instead of a vortex-antivortex pair this pro- 
cedure now creates one long Abrikosov flux line through 
the wire's bulk. The top of the energy barrier is reached 
when the flux line is along a cross-sectional diameter of 
the wire, see Fig. 2. The energy of this critical flux 
line is Eq cx Ly. The rate of change in Ux via ther- 
mal fluctuations is proportional to the Boltzmann factor 
exp(— i?o/T) and thus decreases exponentially with Ly. 
At zero temperature, when spontaneous topological tran- 
sitions have to be through tunneling, the suppression is 
even stronger: a tunneling path is now a worldsheet in 
the Euclidean spacetime, and the tunneling rate goes as 
an exponential of L^. 




FIG. 2. Critical fluctuation in a solid type-II superconduct- 
ing torus. Shown is a cross section of the torus. The critical 
fluctuation in this case is an Abrikosov flux line lying along a 
cross-sectional diameter. 



What we need for the case of a film is a definition of 
the rate of topological transitions that would make no 
mention of a critical fluctuation. This requirement is in 



fact familiar from studies of topological transitions in the 
electroweak theory, where depending on the temperature 
one may or may not have a critical fluctuation to expand 
about. A general definition of the rate in that case is ob- 
tained by considering topological transitions as diffusion 
(or random walk) of the winding number |^ . The rate of 
the transitions is simply the diffusion rate. Here we con- 
struct a similar definition for toroidal superconducting 
film. 

As we already mentioned, to calculate the rate for the 
film it is sufficient to consider the "mathematical" torus, 
on which the classical vacua are given by (||)-(0) and are 
exactly degenerate. To describe diffusion of the winding 
numbers, we need to generalize their definition so that 
it will apply away from the vacuum configurations. This 
generalization is not unique, but the result for the rate 
will be the same as long as the newly defined winding 
numbers are equal to Ux and Uy on the classical vacua 
(^)-(0). A suitable definition is 

a. = ^Jd^.Ax. (14) 

Note that the winding numbers ax and ay are noninteger 
away from the classical vacua. 

Diffusion of ax, ay is due to diffusive motion of vor- 
tices. We assume that the sample is homogeneous enough 
so that most of the vortices are not pinned. Translational 
motion of vortices is semiclassical, so we can define the 
diffusion rates from the classical equilibrium correlator 

{[ax{t)-ax{Q)f)^2Txt (16) 

and a similar one for ay. The linear dependence on time 
on the right-hand side is characteristic of diffusion (in 
the absence of external forces), and F^; is the definition 
of the rate. Eq. ([l^ ) applies at times large compared 
to some microscopic time characterizing interactions of 
vortices with the heat bath. 

The precise meaning of the classical averaging in ( |l^ ) 
is as follows. For each set of initial conditions (for the full 
fields A and we compute ax{Q), then evolve the sys- 
tem until time t, and compute axit). The square of the 
difference is then averaged over all initial conditions, us- 
ing the Boltzmann distribution for those. At this point, 
we should remember however that the system (^) is not 
isolated but evolves under the influence of a heat bath. 
The heat bath is comprised by all degrees of freedom 
not explicitly present in @j — specifically, those associ- 
ated with electrons. So, the requisite evolution equation 
includes a random (Langevin) force, and we need to av- 
erage over realizations of that force as well. 

Before we proceed, it is convenient to recast the def- 
inition of the rate into a different form, which is more 
convenient for actual calculation. The procedure is com- 
pletely standard. First, the left-hand side of ( p^ ) is triv- 
ially rewritten as 
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dt' / dt" {a^{t')a^{t")) . 

"'0 



(17) 



The correlator of time derivatives in (17) is an equilib- 
rium correlator and thus depends only on the difference 
t' — t" . We assume that the corresponding correlation 
time is finite (this assumption can be verified in our spe- 
cific case). Then, at large t the integral ( p7| ) is well ap- 
proximated by 



(18) 



which allows us to rewrite the definition ( p^ ) of the rate 
Vx as 



= 2 / dT{ax{T)ax{0)) 



(19) 



As we will now show, the rate can be found explicitly 
by a simple calculation based on the picture of diffusing 
vortices. 



IV. CALCULATION OF THE RATE 



When a vortex crosses line y — b, the line integral 



Cb = Axdx 

Jy=b 

changes by the amount of the vortex flux, i.e. 
Cb^Cb± 27r/g , 



(20) 



(21) 



the sign depending on which direction the vortex is 
headed. If a vortex moves the entire length Ly (in the y 
direction), it crosses all such lines and changes ax, which 
is essentially the average of gCb/'^T^ over &, by ±1. So, 
if a vortex moves a distance Ay, it changes Ux by the 
amount 



Ay 



(22) 



Taking into account all the vortices (of which there are 
Ny) and antivortices (of which there are Na), we then 
obtain the time derivative of ax as follows 



y \v=l a=l / 



(23) 



We now substitute this expression into the formula ( |l9| ) 
for the rate and assume that, because the vortices are 
well separated, the velocities of different vortices are un- 
correlated. We obtain 



2LI 



dT{y{T)y{0)) 



(24) 



The correlator of velocities in (^J) is computed using 
the equation of motion for a single vortex. We use a 
simple Langevin equation of the form 



Mr + rir = f{t) , 



(25) 



where M is the mass of a vortex, rj in the viscosity coeffi- 
cient, and f{t) is a random force, which we assume to be 
Gaussian white noise; r is the position vector of the vor- 
tex, r — {x,y). The condition of applicability of ( p5| ) is 
that the response of the electronic subsystem to changes 
in ip and A is local; otherwise, there would be a non- 
local response kernel instead of the single coefficient 77. 
The response is local when the mean-free path Ztr of the 
electrons is much smaller than the characteristic length 
scale from which 77 receives the main contribution. As we 
will see in Appendix, the latter length scale is the coher- 
ence length of the superconductor f , so the condition of 
applicability of (|25|) is 



^tr <C 



(26) 



i.e. the superconductor should be sufficiently "dirty" . 

Calculation of rj had a long history and has eventually 
been achieved on the basis of microscopic theory pO[ . 
It is more or less straightforward, though, to obtain an 
estimate, so we present it here. (We assume that the con- 
dition ( P6[ ) is satisfied.) A moving vortex will constantly 
transfer parts of its kinetic energy to the electrons, which 
they will dissipate in collisions with lattice impurities. 
There are two mechanisms of dissipation m^. One is 
Joule heat, which dissipates an amount aE^oi energy 
per unit time per unit volume; here a is the normal con- 
ductivity of the metal, and E is the electric field created 
by the vortex motion. The other mechanism is related to 
response of the electrons to changes in the magnitude of 
-0; it dissipates an amount of order a((9t'0o)^Ttr, where rtr 
is the electronic mean-free time, and a is the parameter 
from (^). These two amounts are typically of the same 
order of magnitude, except at temperatures close to crit- 
ical, where the second amount is small. We estimate 
created by a moving vortex in Appendix. This allows us 
to estimate rj from 



rjv 



a / d^xE^ 



(27) 



where v is the vortex speed. The vortex mass M can be 
estimated from 



1 



1 

8^ 



d-^xE^ 



(28) 



In Appendix, we find that the integrals in (p7|)-(28) are 
saturated at distances r ~ ^ from the vortex center. Cu- 
riously, in our final formula for the transition rate, 77 and 
M will appear only via the ratio 



7 = rj/M ^ a . 



(29) 
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Note that this ratio grows with a, i.e. it is larger in a 
purer metal (which is still "dirty" , though, in the sense 
of (^6|)). Physically, this is because electrons in a purer 
metal more readily accept energy from a moving vortex. 
From (pa), it follows that 



(mm) - (y')exp(-7|r|) 



(30) 



where 7 = r]/M, and {y'^) can be determined by equipar- 
tition: 



Assembling the pieces together, we obtain 
„ TN, + Na 



LI 



(31) 



(32) 



A striking feature of this result is that it does not con- 
tain any volume enhancement: although there is a macro- 
scopic factor of [N^ -\- Na), it is essentially canceled out 
by the inverse powers of Ly. The physical reason be- 
hind this suppression is the extremely long time it takes 
a vortex to circumnavigate the torus: diffusion through 
a distance of order Ly requires time of order Ly . 

The total number of vortices and antivortices is deter- 
mined by the Boltzmann distribution: 



2V 



exp[-/3(Fo + p^/2M)]d?p 



-e-^«/^Mr 



(33) 



where V = LxLy is the total 2d volume and Fq is the 
free energy required to create a vortex. Using _Fo instead 
of the vortex energy takes into account thermal popula- 
tion of the vortex's internal states. Substituting ( |33[ ) into 
(p3), we finally obtain 



l^ 

m-j Ly 



.-Fo/T 



(34) 



This is the rate of transitions that change ax- The rate 
of those that change Uy is obtained by interchanging L^ 
and L„. 



V. DISCUSSION 

As we have already mentioned, for superconducting 
film the exponential factor in ( ^^ can be made practi- 
cally as small as one wishes, because Fq grows linearly 
with the film's thickness. So, a thick film on the surface 
of a torus or, as the limiting case, a solid superconducting 
torus such as shown in Fig. 2 provide quantum memory 
that is stable against thermal fluctuations. We propose 
the following way to write to and read from this quantum 
memory. 



Because magnetic field trapped in the hole of a su- 
perconducting torus (or of any other shape with a non- 
contractible loop) carries energy, the torus behaves as a 
giant "atom", in the sense that it has a discrete energy 
spectrum, with different levels corresponding to differ- 
ent values of Ux- We can write the absolute value of the 
energy difference between levels with Ux = ni > and 
nx = n2 > Til as 



nl), 



(35) 



where R is of order of the linear size of the system (cf. 
Sect. 2) and may depend (presumably weakly) on ni 
and ri2. (We have restored h in this formula.) The cor- 
responding electromagnetic wavelength is 



A = 



27raEM^? 



(36) 



where a-EM is the fine-structure constant. For R of or- 
der of a few cm, and ni^2 ~ 1, the wavelength given by 
(^ ) is in the millimeter range. It is possible that one 
will be able to write a linear superposition of quantum 
states to this device by subjecting it to a pulse of radia- 
tion of frequency w in a resonant cavity, similarly to how 
one induces Rabi precession in atoms. One may be able 
to read from this quantum memory by transferring the 
linear superposition to radiation field in a high-Q cav- 
ity, as was done for atoms in the experiment of ref. p]| . 
Unlike a single photon in a cavity or an excited state of 
an atom, the basis states in our case are macroscopically 
entangled, so this device will be able to store the linear 
superposition for a much longer time. 

If one wants to operate the read and write cavities at 
their principal resonant frequencies and use single-photon 
transitions, at least one of the dimensions of each cavity 
should be of order A. We propose to use, as quantum 
memory, a loop of superconducting wire, such that the 
cross-sectional diameter of the wire is of order A, while 
the size of the loop itself is large enough for R to be 
on the order of centimeters. Only short arcs of the loop 
need to pass through the write and read cavities. This 
arrangement corresponds to Ly ^ X in our formulas. We 
expect that the effective size of the interaction region, for 
interaction between cavity photons and the wire, is also 
of order A, and hence of the same order as Ly. In this 
sense, the interaction is nonlocal, so writing time may 
be not exponentially large. At the same time, Ly is still 
macroscopic, so the rate of thermal transitions changing 
Ux is suppressed. 

We leave calculation of the rate of topological tran- 
sitions induced by a radiation field for future work and 
turn, briefly, to systems with non-Abelian anyons. The- 
oretically, a diverse set of manipulations on degenerate 
states is available for some of these systems It has 
been argued that non-Abelian anyons are realizable as 
excitations of the Pfaffian state Q . The latter is a quan- 
tum Hall state with a certain type of pairing correla- 
tion between electrons and is closely related to the state 
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proposed in P4 ] as a possible explanation of the experi- 
mentally observed u = "0/2 Hall plateau. With non- 
Abelian anyons, nontrivial topology is not required for a 
sample to have degenerate ground states. It is sufficient 
to "puncture" the surface of the sample with a few local- 
ized excitations (vortices). If the typical distance L be- 
tween these vortices is macroscopic, one expects that the 
zero temperature tunneling between the ground states is 
suppressed exponentially with L 

At finite temperature, in addition to those carefully 
planted vortices there will be a sea of thermally excited 
ones. What will be the rate at which quantum memory 
deteriorates in this case? The exponential Boltzmann 



factor, like the one in (34), should still be present in the 
rate. In quantum Hall samples -Fb will be the larger of 
the free energy required to create a vortex and the free 
energy required to unpin it from lattice defects. Because 
these systems are intrinsically two-dimensional, one can- 
not increase Fq at will. As an estimate of fo, we can 
use the value of temperature corresponding to the onset 
of strong temperature dependence of diagonal resistivity. 
This value can be determined experimentally. According 
to ref. |^3|, it is 100 mK for the v — ^>/2 state described 
in that paper. 

The preexponential factor (prefactor) in the rate will 
be determined by motion of thermally excited vortices 
around a localized one. By analogy with the results of 
the present paper, we expect that a thermal vortex that 
is initially at distance p from the localized one will con- 
tribute an amount of order to the prefactor. Then, 
the prefactor will be proportional to 



pdp/ ~ InL' 



(37) 



where L' is either the linear size of the sample or the 
distance between the localized vortices, so the volume 
enhancement of the rate will be at most logarithmic. 

While this paper was being completed, we have learned 
about a recent proposal to use, as a basis for quan- 
tum computation, current-carrying states in supercon- 
ducting loops with Josephson junctions. The authors of 
ref. [|l5[ propose to obtain linear superpositions of these 
basis states by modulating magnetic fluxes through the 
loops with pulses of external current. This technique may 
work also for the quantum memory device proposed here, 
i.e. one may be able to use an external current instead 
of a resonating cavity to change Ux- We plan to return 
to analysis of this possibility elsewhere. 
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APPENDIX: ELECTRIC FIELD OF A MOVING 
VORTEX 

Electric field produced by a moving vortex determines 
the vortex mass M and the viscosity coefficient 77. Here 
we will compute the electric field produced at large dis- 
tances from the vortex core. We will learn in the process 
that the region away from the core is not where most 
of the energy associated with the electric field is concen- 
trated. This precludes us from actually calculating the 
vortex mass, but we will obtain an order of magnitude 
estimate. 

We begin with a collection of formulas describing a 
static vortex, in notation close to that of ref. . The 
magnetic field of the vortex is in the z direction. We 
consider the extreme type-II case when the penetration 
depth 5 of magnetic field is much larger than the coher- 
ence length ^. For the GL Hamiltonian (0), 



1/5^ = 2g\i^l 
lie = 2a/C . 



(Al) 
(A2) 



When distance r from the center of the vortex is much 
larger than ^, the magnetic field of a static vortex located 
at the origin is approximately 



B{x,y)^^Ko[r/5) , 



(A3) 



where r = (x^+y^)^/^, and Kq is the Macdonald function 
of the zcroth order. At r ^ 0, this magnetic field satisfies 



J^V^B -5 = 



(A4) 



We also recall that at small values of its argument A'o is 
logarithmic: Kq{z) = — Inz-I-O(l). 

Now suppose the vortex moves through the origin with 
velocity v, which lies in the x-y plane. The rate of change 
of the magnetic field is 



^tB{x,y) = -^^■VB{x,y) 



(A5) 



The changing magnetic field produces an electric field, 
which is related to dtB via one of Maxwell's equations, 
(V X E)z — —dtB/c. The general solution to this equa- 
tion in our case is 



E(a;, y) = -c-\w x e,)B{x, y) + V/(a:, y) 



(A6) 



where / is so far an arbitrary function. We fix / from the 
condition that V-E = 0. This condition expresses the ab- 
sence of charge separation inside the material; we expect 
it to hold to a good accuracy because charge separation 
in a metal is ass ociated with a large (plasmon) frequency 
gap. Using (A4), we then obtain, at large distances from 
the core, 

E(a;,y) = -c"^(v x e,)B{x,y) 

+ c-^S^{vydx - Vxdy)VB{x,y) . (A7) 
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At small r, the first term here goes as Inr, but the second 
term goes as l/r^. When the second term dominates, 

Kinetic energy of the vortex is 



K = -Mv^ 
2 



1 

8^ 



(A8) 



For the field ( A_7), the integral in (|A8D d iverges at small r, 
due to the singular second term in ( [AT] ). This means that 
the main contribution to the mass comes from the core of 



the vortex, where ( |A7D does not apply. Nevertheless, we 
can obtai n an order of magnitude estimate for the mass 
by using (A7) and cutting of the divergence at distances 
of the order of the core radius, r ^ This gives 



M 



(A9) 



where d is the thickness of the film. The second estimate 
in ( |A9| ) uses the upper critical field Hc2 ~ $o/^^, where 
$0 = TTc/e is the flux quantum. From eq. (p9[), we can 
now obtain an estimate for the viscosity coefficient rj: 



1 
d 



ec 



(AlO) 



which is in agreement with the results of calculations 
based on microscopic theory |ic| ]. 
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